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Abstract 



The mapping of the Wigner distribution function (WDF) for a given bound- 
state onto a semiclassical distribution function (SDF) satisfying the Liouville 
equation introduced previously by us is applied to the ground state of the 
Morse oscillator. Here we give results showing that the SDF gets closer to the 
corresponding WDF as the number of levels of the Morse oscillator increases. 
We find that for a Morse oscillator with one level only, the agreement between 
the WDF and the mapped SDF is very poor but for a Morse oscillator of ten 

levels it becomes satisfactory. 
PACS numbers: 03.65.-w; 03.65.Sq 
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I. INTRODUCTION 



The Wigner description in phase-space [1] provides a tool for comparing quantum and 
classical dynamics [2]. In a previous paper [3] we introduced a mapping relating the 
Wigner distribution function (WDF) corresponding to a given wave function, solution of 
the Schrodinger equation to a semiclassical distribution function (SDF) satisfying the clas- 
sical Liouville equation with the same potential. So far this mapping was applied to the 
ground state of the square well [3], the Poschl- Teller potential [4] and to a modified harmonic 
oscillator potential [4]. 

Here this mapping is extended to the ground state of the unidimensional Morse oscillator. 
We vary the depth of the potential in order to study the effect of the level density on the 
mapping. 

The Morse oscillator [5] plays an important role in many areas of Physics, because it is 
a good approximation to diatomic molecular potentials [6]. Today the Morse potential is 
used in molecular spectroscopy, even for polyatomic molecules [7]. Some authors have used 
the Morse potential for studying the semiclassical limit of Quantum Mechanics [8]. 

The unidimensional Morse oscillator in phase-space has been already studied by many 
authors. The classical motion [9] as well as the quantum representation [6] are well known 
for this potential. 

In Sec. II of this work we summarize the method developed in Ref. [3] for mapping WDF 
into SDF. In Sec. Ill we study the Morse oscillator and in Sec. IV we present our results 
and conclusions. 

II. DESCRIPTION OF THE METHOD 

In this section we present a short derivation of the method used for the mapping [3]. We 
start with the quantum Liouville equation for the WDF p: 

^ + ^Yq + J rK ^P-P)P^P'^)dP ' = 0, (1) 

3 



where the kernel K is given by 



(2) 



V(q) being the potential. The corresponding classical Liouville equation may be written 
similarly 



dp c p dp c 



dt m dq 



J K c (q, p - p')p c (q, p', t)dp' = 0, 



(3) 



where 



i f dv r i 1 dV dV 

K ^> p - p ' ) = -hJM 6XP [h iP - P> \ V = -^d-V 5{P - P) 



dq dq dp 

We relate solutions p and p c of Eqs. (1) and (3) through the integral equation 

/oo r 
dt' / dq'dp'G c (q,p; q',p', t - t') 
-OO J 

x / dp"[K(q',p'-p") - K c (q',p'-p")]p(q',p",t'), 



(4) 



(5) 



where G c is the retarded Green's function corresponding to the classical Liouville equation 
(3). In fact Eq. (5) is the defining equation of p c since the WDF p is already determined by 
fixing the wave function solution of the Schrodinger equation. Using the differential equation 
satisfied by G c [3] one may easily verify that p c satisfies Eq. (3) provided p obeys Eq. (1). 

As K c is the term of lowest order in the expansion of if in a power series of h, p — p c 
contains only corrections of the first and higher order terms in h. 

It can be shown from results given in Ref. [3], for a WDF generated from solutions of the 
Schrodinger equation for which the initial wave function at the time t — to is given; assuming 
p(q,p,t) to vanish for t < t , that p c from Eq. (5) gives simply the classical evolution of 
p(q,p,t ), that is 



p c (q,P,t) = p(Q(q,p,t - t),P(q,p,t -t),t ), t > t . 



(6) 



Here (Q(q,p,T), P(q,p,r)) describes the classical trajectory in phase-space of a particle 
subject to the potential V(q) which at the time r = occupies the point (q,p), that is 



q = Q(q,p,0), p = P(q,p,0). (7) 

In what follows we shall consider mostly the mapping of the time independent WDF 
p(q,p) corresponding to a bound state. In this case, from Eq. (5) one obtains [3] the 
following prescription for p c , 

p c (q,p) = lim e f dr e er p(Q(g,p, r), P(q,p,r)), (8) 

which is also time independent. Here e €T is a convergence factor introduced explicitly in the 
Green's function. 

It can be shown [3] that for points (q,p) on closed classical orbits the expression given 
by Eq. (8) is equivalent to 

1 r° 

p c (q,p) = 7^7 r/ dtp(Q(q,p,t),P(q,p,t)), (9) 

T(q,p) J-T(q,p) 

where T = T(q,p) is the period of the orbit. For open trajectories we have T — > oo and as 
for bound states the integral in Eq. (9) converges, we get p c = 0. As it will be discussed 
further, in the calculation of averages of physical quantities the quantity of interest is rather 
T(q,p)p c (q,p) which does not vanish for the open trajectories. 

We observe here that one may verify directly that the time-independent Liouville equa- 
tion is obeyed by applying the operator pd/dq — (dV / dq)d / dp on the right hand side of 
Eq. (8) or (9) and using the fact that p(Q(q,p, —i),P(q,p, —t)) satisfies the classical Liou- 
ville equation (3). 

An important property of the stationary SDF is that it is constant along the classical 
trajectories [3] which can be seen by noticing that for any two points (q,p) and (q',p') on 
the same path in phase space one has a time interval A such that 

Q(q', p', t) = Q(q,p,t + A), P(q', p', t) = P(q, p,t + A). (10) 

(A is the time which the particle takes to go from the point (q,p) to the point (q',p') in 
phase space or vice versa). Thus we get 
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Pc(q',p') = \ J° T dtp(Q(q', P ',t),P(q', P ',t)) 

= ^ £ T dtp(Q(q,p,t + A),P(q,p,t + A)) 

= pc(q,p), (11) 

where we used the periodicity of the integrand in the last step. 

Let (Q(E,i/,t),P(E,v,t)) represent, for t in the interval -T(E,v)/2 < t < T(E,u)/2 
(T may be infinity), the points of a trajectory in phase-space corresponding to energy E and 
fixed value of v. Here the discrete parameter v {y — 1, • • • , z/ max ) distinguishes between the 
different trajectories corresponding to the same value of E. By varying E and v one covers 
the entire allowed phase space. Thus we may consider the transformation of points (E, u,t) 
in the appropriate domains D u of the energy-time space onto the points (p, q) of the phase 
space 

q = Q{E,v,t), p = P(E,v,t). (12) 
Let us define the function 



R, 



{E,v)= / dtp(Q(E,v,t),P(E,v,t)). (13) 

J-T(Em)/2 



-T{E,v)/2 

According to Eq. (9) one has 

R c {E(q,p),u) 

P < ( "- P) = T(E( q , pM - (14) 

by choosing the value of v appropriate to the trajectory containing the point (q,p) as follows 
from Eq. (12). We observe here that the relationship (13) between p(q,p) and R c {E,v) is 
analogous to that between R(q), the probability density in coordinate space and the WDF 
p(q,p) [10]: 

R(q) = Jp(q,p)d P . (15) 

The average of any Weyl function [10] 0(q,p) corresponding to a certain operator O 
may be written 
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(0) = j0(q,p)p(q,p)dpdq = J2j v o(E, v, r)r(E, u, r)dEdr, 

where we introduced the functions 

o(E,u,t) = 0(Q(E,u,t),P(E,u,t)), 
r(E,u,r) = p(Q(E,u,t),P(E,u,t)), 

and we used the fact that the Jacobian J of the transformation (12) is unity [3] 

dP dE dQ dE 



(16) 

(17) 
(18) 



J 



dP dQ 
dE dE 



dP dQ 
dt dt 



dE dp + dE dq 



(19) 



In the special case in which the Weyl function O is a constant of motion depending on 
(q,p) only through the energy E(q,p), we obtain from (16), (18) and (13) 

(0)=J2 [ o(E, u)r(E, v, r)dEdr = V / o(E, v)R c {E, v)dE. (20) 

In particular the normalization condition for R C (E, v) 

I R c (E,u)dE = l, (21) 

v J 

follows from Eq. (20) by taking 0(q,p) = 1 since the Wigner function p(q,p) is normalized. 

As p c is an approximation correct in zeroth order of the expansion in powers of h of the 
Wigner function p, the average 



(0)c 



0(q,p)p c {q,p)dqdp, 



(22) 



is also correct in the same order in H. The expression (22) may be also written, by making 
use of the transformation (12), as 

(0> C = W o(E,u,t) R r jf u) dEdt = J2 [ o(E,u)R c (E,u)dE, (23) 

where we used Eqs. (14) and (17) and d(E, v) denotes the average 

o(E,v) 



T(E,u)/2 dt 

o(E,is, t)———. 



(24) 



Thus the average {0) c is obtained by replacing o(E, u, t) in Eq. (16) by the average (24). 
If, for fixed E, o(E,u,t) depends weakly on t, the average (0) c is expected to be a good 
approximation. 



III. WIGNER DISTRIBUTION FUNCTION FOR THE MORSE POTENTIAL 



In 1929 Morse [5] suggested the potential U (l — e~ ar ) 2 for studying diatomic molecules. 
The Schrodinger equation for this potential does not have an exact solution, but for the one 
dimensional case the problem can be solved analytically [11,12]. 

In order to obtain the Wigner Distribution Function (WDF) [1] 



/oo 1 I 

dve^(q--v,t)r(q + ^v,t), 



(25) 



where ip is the solution of the Schrodinger equation, we need the eigenfunctions for the 
one-dimensional Morse potential 



V(x) = D(l - e' ax Y, 
a and D being constant parameters. Starting with the Schrodinger equation 

h 2 d 2 



2m dx 2 



+ D(l-e 



—ax\2 



we introduce the dimensionless parameter A 

V2mD 



A 



ah 



and the dimensionless coordinate 



q = ax, 



obtaining an eigenvalue equation depending only on one parameter 



1 d 2 



where 



£ 

~D 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



This equation is solved most conveniently using the variable 
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£ = 2Ae 9 , -oo < q < oo. 



(32) 



The eigenfunctions and eigenvalues are [6] 



VUO=^A)e-^ A -"-^- n "(0, 



1 / 1 
1 -2A( n+ 2 



(33) 
(34) 



where the quantum number n takes the values 



n = 0,l,2,...[A--]. 



(35) 



Here [A — 1/2] denotes the largest integer smaller than A — 1/2, is the polynomial [6] 

(36) 



3=0 



n ( n + 2s\ 



n-j 



J- 



and the normalization factor N is given by 

(2A-2n- l)r(ra+l) 



N(\,n) 



(37) 



r(2A-n) 

where the normalization condition / i/j*ip n dq = 1 is assumed. In Appendix A, Eqs. (33) and 
(34) are derived. 

From Eq. (34) one verifies that, for A 3> 1 and n<C A, the energy spectrum of the Morse 
oscillator is written approximately S n ~ hcu (n + 1/2), which is the spectrum of a harmonic 
oscillator with frequency 

2D [2D x + a 2 
u = — = a\ — = An — . 
n\ V m m 



(38) 



As the semiclassical distribution functions p c for the harmonic oscillator coincide with 
the Wigner functions p we expect that if A ^> 1 p c does get close to p for the low lying levels. 
Thus it may be appropriate to use, instead of the variable q and its canonical conjugate 
momentum p, the variables which treat harmonic oscillators on the same footing, namely [6] 



Q 



nr. 



X 



-2mD- 
_a?h 2 \_ 



y/Xq, 



(39a) 
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p = 7f- 



(39b) 



The coordinate Q and the momentum P have been used in the figures of Sec. IV. We 
define the dimensionless potential V(Q) 

V(Q) = (hu )- l V{x) = ± (l - e-%y , (39c) 

which is also used in the figures. For the energy levels there we use similarly 

E n = (huo^Sn = (n + - (n + £) . (39d) 

In the special case of our interest, namely the ground state, n = 0, the wave function is 
given by 

"2A- 1' 



MO 



r(2A) 



*\ 1 « 

? 1 e 2, 



(40) 



In order to calculate the Wigner function replace ip by ^ n (2Ae 9 ) into Eq. (25), obtaining 

/oo 
dvM^e^m^Xe^) ei*°. (41) 

Introducing the new integration variable 



Eq. (41) becomes 



t = e 2 , 



P ( n X \q,p) = (Kh)- 1 f-1 ^(er)r 2 ^ 1 -, 



(42) 



/o \ T J r 
where £ is given by Eq. (32). Substituting ip n in Eq. (43) by the expression (33) p, 
becomes [6] 



(43) 



(A) 

v. 



where 



-\ -i 



[N(X, n)] 2 ^ 2 - 1 £ E K\ n, r)b(X, n, s)C + > 'K a _ r+MpR - 1 (0, (44) 

r=0 s=0 



b(X,n,j) 



^ 2X-n-l \ (-1) 
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(-1V 



(45) 



and K v (£) is denned by [6] 

dr 



2 Jo T 



v being a complex variable. In the particular case in which n = we get 



P £\ q , P ) = (^-y 1 ^e^K^AZ). (47) 
The numerical method we used to calculate the function K v (£) will be given in Appendix 

B. 

The symmetries obeyed by the function K v are 

^(0 = \ f 'e-^+^rTv = (^(0)*. (48) 
2 jo r 

= £ r e-fr + &T- v - = (49) 

2 JO T 

where the last step is obtained by making r — > r _1 . From Eqs.(49) and (48) we get 



K^ a+ib = K* a+ib . (50) 



Thus Eq. (44) may be written 



P ( *\q,p) 



1 [iV(A,n)] 2 e 2A - 2n - 1 EEKA,n,r)6(A,n, S )r + ^e (^+^(0) , (51) 

/ J r=0s=0 

showing explicitly that p$ is real. 

In order to calculate the semiclassical distribution function p c (q,p) given by Eq. (9) or 
(14), we need the solution of the classical equation of motion 

mx = -D-^-il -e~ ax )\ (52) 
ax 

which has already been obtained exactly [9]. Introducing the coordinate q = ax and the 
variable 9 given by 

9 = u t, (53) 
where u is given by Eq. (38), Eq. (52) becomes 
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The solution of Eq. (54) for e = E/D < 1, where E is the energy associated with the 
trajectory, is 

1 



q(E,t) = \n 



1-6 



1 + \fk sin yl — e{uJot — 9 ) 



The canonical momentum associated with q is 



mq 

V = -5", 

cr 



which using Eqs. (55) and (38) gives 

p(E, t) = h\ ^^costv^eKt-go)] 
1 + v/esinfvl - e(cj ^ - #o)J 

The period associated with the orbit is 

2tt 



T(E) 



(55) 



(56) 



(57) 



(58) 



Analising Eq. (55) one obtains that foiE<^D the trajectories are close to those of a 
harmonic oscillator of frequency uoq. 



IV. RESULTS 

In this section we present the results of our calculations. We calculated the WDF p and 
corresponding SDF p c for the ground state of the Morse oscillator choosing for the parameter 
A the values 1,2,4 and 10 corresponding to oscillators with 1,2,4 and 10 levels respectively. In 
the figures we used throughout the dimensionless coordinates Q and the conjugate momenta 
P (in units of H) defined by Eqs. (39a) and (39b). In Figs. 1,2,3 and 4 the potentials V(Q) 
defined by Eq. (39c) are displayed for A = 1,2,4 and 10 and the energy levels marked. We 
observe here that the potential given by Eq. (39c) is independent of A for Q 2 <^ A. 

Figs. 5,6,7 and 8 reproduce our calculations of the WDF for A = 1, 2, 4 and 10 respectively 
through curves of constant density p(Q, P). The value of p varies from p ps 0.3 to p = in 
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the region of phase-space corresponding approximately to the region of the bound classical 
particles. Then there occur adjacent regions on which p alternates from negative to positive 
values and its magnitude decreases as the region gets farther away from the origin. 

In the case A = 1 the minimum value of p is (in units of hT 1 ) of the order of — 1CT 2 , 
as it can be seen from Fig. 5. This minimum approaches zero as A increases, which can be 
observed from Fig. 7 for A = 4 where the minimum of p is about — 10 -4 . 

The maximum of p, as it can be verified from Figs. 5 — 8, moves from the point (Q,P) = 
(1.2,0) to (Q,P) = (0.3,0) as A increases from 1 to 10 and its value increases slightly as A 
increases. Thus the curve for p = 0.3, present for A > 2, does not occur for A = 1. 

Another feature of the WDF is that the curves of constant density p become more 
symmetric with respect to an axis parallel to the P— axis as A increases, becoming close to 
the form of an ellipse. This tallies with the fact that for Q <C y/X, V(Q) is the potential of 
a harmonic oscillator. 

In Figs. 10,11,12 and 13 we present SDF curves for fixed p c superposed on WDF curves 
with p = p c for comparison. It will be seen that, as A increases the SDF approximation 
improves, which means also that the WDF curves of constant p become closer to classical 
trajectories. In Fig. 9 we plot curves of constant p c for the case A = 1. For this oscillator 
our semiclassical approximation is anomalous, as the value of p c increases from p c = 0.145 
to p c = 0.179 as the classical energy E varies from E = to E = 0.2Qhuj but decreases as 
E increases further. For the other oscillators p c decreases as the energy E increases until 
reaching the value p c = 0. As a consequence for A = 1 one has two closed curves with the 
same p c for 0.145 < p c < 0.179 whereas for A > 2, for each p c one has only one such curve. 
Also for the WDF there is only one closed curve for each p from the maximum value of p 
up to the curve p = 0. 

In Fig. 10 we compare the WDF curves of constant p with the SDF curves for which 
p = p c in the case A = 1. One notices that the discrepancies between both curves are very 
large. In Figs. 11,12 and 13 we make the same comparison for A = 2,4 and 10 respectively. 

One finds that for p < 0.05 both curves are quite similar but displaced from each other. 
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This displacement becomes less pronounced as p gets closer to 0.05 so that the best agreement 
between p and p c is reached for p ps 0.05. Also as A increases the displacement between 
both curves decreases, as it can be seen comparing the oscillators A = 4 and A = 10. For 
p ~ 0.05 the curves of fixed p c compared with the curves of fixed p contain a certain amount 
of distortion which becomes less pronounced as A increases. 

Finally we observe that curves with p c = 0.3 are absent, the largest values of p c being 
0.227, 0.271 and 0.299 respectively for A = 2,4 and 10, values which are reached at the 
origin of phase-space. 

APPENDIX A: EIGENFUNCTIONS AND EIGENVALUES OF THE MORSE 

POTENTIAL 



In this Appendix we solve briefly the eigenvalue equation for the Morse potential (Eq. (30) 
of Sec. Ill) 

d 2 



__ _ a 2 (1 - e~ z f + X 2 e n 



Making the substitution y = 2Ae z , Eq. (Al) is written 



y 



+^-A 2 (l-^) 2 + A 2 e r 



dy 2 dy 



i> n {y) = o, 



(Ai) 



(A2) 



or, 



d 2 Id e n - 1 A 2 A 1 
+ -— + V ' +- - - 

y 4 



dy 2 ydy y 
Assuming a solution of the form 



Mv) = o. 



(A3) 



ipn{y) = e *y s u n {y), 



(A4) 



and writing s 2 = (e„ — 1)A 2 , from Eq. (A3) we get 



yv!' n + (2s + 1 - y)u' n + (A - s - -)u n = 0, 



(A5) 
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where u' n = du n /dy and v!' n = d 2 u n /dy 2 . Making u n = J2m=o a mU m , Eq. (A5) may be written 

as 

n ( 1\ 

[(m + l)(m + 2s + l)a m+ i + A - s - m - - a m ]y m = 0, (A6) 

m=0 V ^ 

which is the differential equation satisfied by the Laguerre generalized functions. Eq. (A6) 
is satisfied if 

(m + l)(m + 2s + 1) 



which for m = is 



2s + 1 

Since for m = n, a m+ i = and a m ^ 0, Eq. (A7) gives 



s — A -\- h 

ai = t 2 - a . (A8) 



2 

Taking into account Eq. (A9) we get 



s = \- n --. (A9) 



2 / IN 1 / 1^ 2 



- f " = -A(,' l+ 2j + ^t' !+ 2j • (A10) 

and the recurrence relation for the coefficients of the power series in Eq. (A7) becomes 

m-n 
(m + 1)(2A — 2n + to) 

Choosing 

^' nlS-I"-!)! ^")- (A12) 



we obtain 



(-ir(2A-n-l)!iV(A,n) 
a - " m \(n _ m ) !( 2A - 2n - m - 1)! ' { ' 



and we get for ip n the expression given by Eq. (33) in Sec. III. 
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APPENDIX B: NUMERICAL EVALUATION OF THE WIGNER FUNCTION 

Here we discuss the numerical calculation of the quantity K v defined in Eq. (46) of 
Sec. Ill, 

A-„ K ) = \ r ,-««*) * (Bl) 

2 JO T 

Here v is complex, 

v = N + 2ik, (B2) 

where N is an integer and k is the dimensionless momentum p/h and according to Eqs.(32) 
and (29) 

£ = 2\e~ q , q = ax. (B3) 

Making the transformation 

r = e u , (B4) 

and considering that only the real part of K v enters into the expression (51) for the Wigner 
function we get from Eq. (Bl), substituting also v according to Eq. (B2), 

1 r°° 

Re[K v {i)} = - e^ coshu e uN cos(2ku)du. (B5) 

2 J—oc 

Let us take initially N = 0, which is the only value needed for the ground state of 
the Morse oscillator. Consider also k ^ as the case k = is calculated separately. For 
convenience we introduce the new variable of integration 

z = 2ku, (B6) 

and use the fact that the integrand is an even function of u, obtaining from Eq. (B5) 

Re [K v (£)] = — / dze-^ osH ^ cosz. (B7) 

The integral in Eq. (B7) is of the form 
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f(z) cos zdz, (B8) 

where f(z) is a positive decreasing function of z. In order to avoid numerical cancellations 
arising from the change of sign of cos(z), we replace this integral by an integral in the interval 
[0, |] of a series of functions. 
We decompose / as follows 

I = h + h, (B9) 

where 

h= 2 f(z) cos zdz, (BIO) 
Jo 

and 

h= f(z)coszdz = - f(y + -)sinydy, (Bll) 
if Jo 2 

where we made the substitution z = y + | in J 2 - The interval of integration of the integral 
J 2 is divided in the set of intervals [2ns, 2n(s + 1)], s — 0, 1, • • • so that J 2 is written 

~ /•2 7 r(s+l) 

/ 2 = -E/ /(?/ + -)sin^. (B12) 

Making now the substitution 

z = y- 2vrs, (B13) 
for the integral in the interval [2irs, 2tc(s + 1)], one obtains for Eq. (B12) 

00 r-2n 1 

J 2 = -E/ /(^ + (2s + -)7r)sin^. (B14) 

However sinz changes sign in the interval [0, 2n] which leads to cancellation errors if / is 
slowly varying in the interval. In order to eliminate the oscillation of sign of the integrand 
we decompose again the integral in the intervals [0, n] and [n, 2n\ obtaining 

J 2 = -E{/ o " [/ (z + (2s + 1) tt) - / (z + (2s + |) tt)] sin zdz} . (B15) 
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As / is assumed to be monotonically decreasing each integrand in Eq. (B15) is now 
always positive. As the domain of integration of the integral Ii is the interval [0, |] we make 
the translation v — z — |, obtaining 

OO „2L 

1 2 = - \ \ \f( v + ( 2s + 1 ) 7r ) - f( v + ( 2s + 2 ) 7r )l cosvdv. (B16) 

Making t> — > — v for the contribution from the interval [— | , 0] in Eq. (B16) and adding 
the contribution from l\ one gets finally 

r°° r% °° 

/ f{x)cosxdx = / dx cos x{ f {x) - J2[f{{2s + 1)ti -x) + f {{2s + 1)ti + x) 
Jo Jo s=0 

- /((2s + 2)tc -x)- /((2s + 2)tc + x)]}. (B17) 

For f{x) monotonic each term of the series contributes with the same sign, however 
errors may arise from the subtraction of the sum of the series from f(x). 

In the general case in which iV > in Eq. (B5), the function f{z) in Eq. (B8) becomes 

f N {z) = ^ cosh ^ cosh , z > 0, N = 0, 1, • • • (B18) 

The function /n(z) has M (M < N) extremes at the real positive points 
z[ N \ z^, ■ ■ ■ , z^\ labeled in the order of increasing magnitude. This function decreases 
monotonically for z > z^ and it may still be useful to apply the decomposition (B17) 
in order to eliminate errors due to cancellations of contributions of opposite sign from the 
integrand fjy{z) cos z. In the general case only the sign of the first sjv terms of the series 
in Eq. (B17) may oscillate, where sn is given roughly by the smallest integer satisfying 
(2sjv + 1)tt > z^\ For s > sn the sign of the terms of the series in Eq. (B17) is always 
positive. 

In fact, one may determine the extremes of /n{ z ) by expressing /at(z) in terms of 
cosh{z/2k) and applying the condition dfjy{z)/dz = 0. By making this substitution one 
gets for Eq. (B18) the expression 

1 [JV/2] _ 

fs{z) = x: i)V 21 ■ y = cosh(-), (Big) 

18 



where [N/2] denotes the largest integer contained in N/2 and 



[N/2] I N ^ 

A i= E 
i=j V 2 V 



Thus we get the extrema of /at (z) as the roots of a polynomial of degree 
For N — 1 the maximum occurs at 



4 X) =2A;cosh- 1 ^ . 



For N = 2 the maximum will be at 



z[ 2) =2A;cosh- 1 Q + + 2 



19 



REFERENCES 

[1] E. Wigner, Phys. Rev. 40, 749 (1932). 

[2] H. W. Lee and M. O. Scully, J. Chem. Phys. 77(9), 4604 (1982). 
[3] G. W. Bund, J. Phys. A28 3709 (1995). 

[4] G. W. Bund and M. C. Tijero, Phys. Rev. A 61, 052114 (2000). 
[5] P. M. Morse, Phys. Rev. 34, 57 (1929). 

[6] J. P. Dahl and M. Springborg, J. Chem. Phys. 88(7), 4535 (1988). 
[7] A. Frank et al, Ann. Phys. (N.Y.) 252, 211 (1996). 

[8] J. P. Dahl, in Semi-classical Descriptions of Atomic and Nucleon Collisions, J. Bang 
and J. de Boer (Eds.) (Elsevier Science Publishers, B.V. 1985). 

[9] W. C. De Marcus, Am. J. Phys. 46, 733 (1978). 

[10] S. R. Groot and L. G. Suttorp, in Foundations of Electrodynamics (Amsterdam: North- 
Holland, 1972). 

[11] M. M. Nieto and L. M. Simmons Jr. Phys. Rev. A 19, 438 (1979). 

[12] D. ter Haar, Problems in Quantum Mechanics, 3 th Edition (Pion, London, 1975). 



20 



FIGURES 
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FIG. 1. The Morse potential V(Q) denned by Eq. (39c) and the corresponding bound-state 
energies E n given by Eq. (39d), for A = 1. 
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FIG. 2. Same as Fig. 1 for A = 2. 
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FIG. 5. Curves of constant Wigner distribution function p(Q, P) in units of % for the ground 
state of the Morse oscillator with A = 1. 




FIG. 6. Same as Fig. 5 for A = 2. 
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FIG. 8. Same as Fig. 5 for A = 10. 
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FIG. 9. Curves of constant semi-classical distribution function p c (Q, P) in units of h^ 1 for the 
ground state of the Morse oscillator with A = 1. For 0.145 < p c < 0.179 one has two curves for 
each value of p c . 



25 



1.5 
1 

P 

0.5 



"-2 2 4 6 

Q 

FIG. 10. Full curves are curves with a constant value of the Wigner distribution function p and 
dashed lines have constant value of the semiclassical distribution function p c for the ground state 
of the Morse oscillator with A = 1. Except for p = 0.2, for each full line with a given p one has one 
or two corresponding dashed lines with p c such that p c = p. 
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FIG. 11. Same as Fig. 10 for A = 2. Except for p = 0.3 and p = 0.25, for each full line one has 
a corresponding dashed line with p c such that p c = p. 
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FIG. 12. Same as Fig. 10 for A = 4. Except for p = 0.3, for each full line one has a corresponding 
dashed line with p c such that p c = p. 




